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SUMMARY 

Computational methods for simulating and predicting progressive fracture 
in fiber composite structures are presented. These methods are Integrated Into 
a computer code of modular form. The modules Include composite mechanics, 
finite element analysis, and fracture criteria. The code Is used to computa- 
tionally simulate progressive fracture In composite laminates with and without 
defects. The simulation tracks the fracture progression In terms of modes 
Initiating fracture, damage growth, and Imminent global (catastrophic) laminate 
fracture. 


INTRODUCTION 

Prediction of progressive fracture of composite laminates Is fundamental 
to developing the methodology for quantifying composite structural durability 
and reliability. NASA Lewis Research Center Is conducting fundamental theo- 
retical and experimental research programs to develop formal methods and pro- 
cedures for determining progressive composite fracture. The theoretical 
studies Include the development of composite mechanics, combined stress fail- 
ure criteria, criteria for Identifying predominant fracture modes and asso- 
ciated fracture surfaces, and the development of Integrated computer codes for 
the computational simulation of progressive fracture In fiber composites., The 
experimental studies Include development of methods for real time ultrasonic 
C-scan of laminates under load In order to detect fracture Initiation, damage 
growth and fracture progression as these events occur. The experimental 
studies also Include development of methods for post-mortem examination of 
fracture surfaces In order to Identify and catalog unique fracture surface 
characteristics associated with dominant fracture modes. 

The theoretical studies previously mentioned led to the development of 
computational methods for assessing composite structural behavior (Integrity, 
durability and damage tolerance). A major portion of the computational methods 
development Includes the computational simulation of progressive fracture In 
composites which Is the subject of this paper. This simulation Is described 
In terms of (1) the various composite scales (as defined later) and types of 
fracture modes that occur In these scales, (2) an Integrated computer code to 
perform the simulation, and (3) application of the code to predict progressive 
laminate fracture. 


FUNDAMENTAL CONSIDERATIONS 

Progressive fracture In fiber composites Is a dynamic (real-time) physical 
process. This process Is an accumulation of multiple events triggered by the 
various stress states and their attendant failure modes present In fiber com- 
posites. The various failure modes take place at different scales (levels) In 


the composite. The different levels Include: (1) sub-micro scale (Intra or 

within each constituent and the Interphase), (2) micro scale (Inter-con- 
stituent), (3) macro scale (Individual plies and Interplles), (4) laminate 
scale (all plies In a line through the thickness of the laminate), (5) local 
scale, a small region defined by neighboring points In the laminate plane (for 
example, nodes connecting a finite element) and, (6) global scale (the struc- 
tural component). 

Because of the several failure modes In each scale, fracture can occur as 
single events at different times In different locations In the composite 
structure. Or It can occur simultaneously as multiple events In one or more 
locations In the composite structures. In view of this. It Is convenient to 
discuss progressive fracture using the following definitions: (1) failure 

Initiation - when the stress field Induces a failure In the sub-micro scale, 

(2) defect or flaw formation - when one or more sub-micro scale failures 
coalesce to form a micro scale failure, (3) ply failure - when several micro 
defects coalesce to form defects In the macro scale (transply cracks), 

(4) delamlnatlon - when a defect Is formed between adjacent plies, (5) laminate 
failure - when macro defects (transply cracks) coalesce to form a through-the- 
thlckness laminate defect or crack, (6) progressive fracture - when laminate 
defects progress In the plane of the laminate to form defects In the local 
scale, and (7) global fracture - when the composite structure fractures In one 
or more parts. Events occurring In Items (1) to (6) are usually combined and 
are collectively called laminate or composite damage. 

Computational simulation of progressive fracture In composites as a real- 
time dynamic event has not yet been Investigated to the author's knowledge. 
Also, computational simulation of progressive fracture as an equivalent static 
process starting from the sub-micro scale has not yet been Investigated. The 
computational simulation presently pursued at Lewis Is based on the macro scale 
and It assumes that progressive fracture Is an equivalent static process. In 
this computational simulation, macro scale defects (transply cracks) are formed 
when the stress field Induces failure In one of the ply failure modes (longi- 
tudinal tension or compression, transverse tension or compression and Intra- 
laminar shear) and/or Induces Interply delamlnatlon. The corresponding ply and 
interply strengths are predicted using composite micromechanics. The composite 
micromechanics account for void and hygrothermal effects on the constituents 
and, therefore, on the ply and Interply strengths. Composite micromechanics 
Integrate, In part, the events occurring at the sub-micro and micro scales. 


COMPUTATIONAL SIMULATION USING COOSTRAN 

Research activities on progressive composite fracture at Lewis during the 
last ten years have culminated In the development of the COOSTRAN computer code 
(ref. 1). COOSTRAN (Composite Durability Structural Analysis) has been speci- 
fically developed for the computational simulation of progressive fracture In 
fiber composites. 

COOSTRAN Is a modular program (fig. 1) that does quantitative calculations 
to predict defect growth and progressive fracture In composite structural com- 
ponents. Capabilities of COOSTRAN Include determining the durability of com- 
posites with and without defects, determining structural responses due to 
mechanical and thermal loads, accurate prediction of stress states near 
defects (stress concentrations), and prediction of ply and laminate level 
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failure and fracture. The modules comprising CODSTRAN are: (1) the executive 

module, containing communication links to all other modules; (2) the I/O 
module; (3) the Analysis module; (4) the Composite Mechanics Module; and 
(5) the Fracture Mechanics module. 

The Analysis module Is NASTRAN (ref. 2) and Is used to calculate both 
near-field and far-fleld stresses In a finite element model of the structural 
component or specimen. The Composite Mechanics module (MFCA) (ref. 3) gen- 
erates laminate properties from constituent properties (composite microme- 
chanics) and uses Intraply failure and Interlaminar delamlnatlon criteria to 
check ply and interply failure, respectively. The Fracture Mechanics module 
Is able to account for both ply level fracture and laminate fracture. The 
modified distortion energy principle and/or a general quadratic surface fit are 
used to indicate combined ply level fracture (ref. 4). Laminate level failure 
Is based on concepts described In references 5 and 6. It Is assumed to occur 
when all the plies In the laminate have failed. 

To computationally simulate progressive fracture, COOSTRAN uses an Iter- 
ative procedure whereby a load Is applied to the finite element mesh of the 
structure being modeled. The response of the structure to the load can be no 
damage, damage, or destruction of an element(s) (local fracture). Based upon 
this response, the load Increment Is updated as follows: (1) If no damage Is 

predicted, the load Is updated by some predetermined load Increment; (2) If 
elements are damaged or destroyed (defect growth or local fracture), the same 
load Is re-applied with reduced material properties assigned to the damaged 
elements. Destroyed elements are purged from the finite element mesh, effec- 
tively defining progressive fracture. This load Is maintained until equilibrium 
In the structure Is achieved. Equilibrium Is defined as the point where the 
structure, with Its updated geometry and modified material properties, can 
sustain the applied load without the occurrence of further damage. This Iter- 
ative procedure is continued through load Increments until global fracture of 
the structure occurs. This Iterative procedure Is depicted schematically In 
figure 2. 


PROGRESSIVE FRACTURE - TYPICAL RESULTS 

CODSTRAN has been used to computationally simulate progressive fracture 
In a composite laminate with a center notch depicted schematically In figure 3 
(ref. 1). These laminates were made from Thornel 300 graphite fiber In an 
epoxy matrix (T300/E). The laminate configuration Is [0/+30/0/-30/0]2 S with 
about 0.005 In. ply thickness. These types of laminates were loaded to frac- 
ture In axial tension. The finite element model used In the simulation Is 
shown In figure 4. 

Progressive damage and fracture results obtained are summarized In figure 
5 at several load levels Indicated as a percentage of the fracture load. First 
damage Indication appeared at 33 percent of the load. The damage continued to 
grow and extend until the load was Increased to 100 percent of the fracture 
load. At this load CODSTRAN looped through several Iterations. The first 
Iteration, caused the extended damage shown In the first 100 percent schematic 
In figure 5. The last Iteration destroyed all elements to the right of the 
cracks, Indicating laminate global fracture with some bending as shown In the 
last schematic In figure 5. 
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CODSTRAN has also been used to study the damage growth and progressive 
fracture of laminates without defects, with center slits and with center holes 
(refs. 7 and 8 ). Progressive fracture results from this study are summarized 
In figure 6 at zero load and at fracture load, first Iteration. As can be seen 
the progressive fracture Is about the same for the laminates with the slit and 
with the hole. It Is Interesting to not** that progressive fracture Initiated 
at the center of the laminate without defects and advanced In generally similar 
directions as that In the laminates with the defects. 

CODSTRAN keeps records of all the modes that Initiate fracture at the 
macroscale (ply and interply) levels. A typical output of these records Is 
summarized In Table I for the laminates shown In figure 6 and for different 
laminate configurations [±e s ]. These results provide considerable details 
with respect to weak modes In the laminate and with respect to Its structural 
Integrity and/or damage tolerance. In addition to records of fracture modes, 
CODSTRAN keeps records of plies destroyed, elements destroyed and nodes which 
do not connect elements. All this Information Is necessary to track the damage 
(defect growth and progressive fracture) at the various scales within which It 
occurs. 

In order to assess the fracture toughness and service life of the laminate, 
CODSTRAN continued to track the crack or defect opening In either the local 
scale (between two adjacent nodes) or global (overall displacements) scale. 
Typical results for local crack opening displacement In laminates with center 
slits are shown In figure 7 for two laminate configurations. The results for 
the [±45] s laminate show unbounded crack opening displacement for the same 
load and, therefore, Imminent global (catastrophic) fracture once progressive 
fracture began. The results for the [ 0)4 laminate showed several Increases 
In load were required prior to Imminent global (catastrophic) failure. 


CONCLUSIONS 

Progressive fracture In fiber composites can be computationally simulated 
using an integrated computer code such as CODSTRAN. The computational simu- 
lation tracks failure Initiation, defect growth and damage, at the various 
scales In which these events occur. The computational simulation of pro- 
gressive fracture provides extensive detailed Information which can be used to 
detect failure Initiation modes, damage growth (magnitude and direction), pro- 
gressive fracture direction, and imminent global fracture. All this Infor- 
mation makes It possible to computationally assess the structural Integrity, 
damage tolerance and service life of fiber composite structures. 
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TABLE I. - FRACTURE MODE* OF [±e] s G/E LAMINATES 
(PREDICTED BY CODSTRAN) 



Ply configuration: 

[±e] s ; 

e In degrees 

Notch type 

0 

3 

5 

10 

15 

30 

45 

60 

75 

90 

Unnotched 

LT 

LT 

LT 

LT 

I 

S 

I 

TT 

TT 

TT 

solid 


S3 

S3 

S3 

S 


S 




Notched 

LT 

LT 

LT 

S 

S 

I 

I 

I 

TT 

TT 

thru slit 

s 1 

S 

S 



S 

s 

TT 











S2 



Notched 

LT 

LT 

LT 

S 

1 

I 

I 

B 

TT. 

TT 

thru hole 

s 

S 

S 


19 

S 

s 

19 




1 

1 




HI 


TT | 

■I 




a LT = Longitudinal tension 
TT * Transverse tension 

S = Intraply shear: 1) Intraply shear occurring 

around notch tip during 
progressive fracture 

2) Minimal Intraply shearing 
during fracture 

3) Some Intraply shear occurring 
near constraints (grips) 

I = Interply delamlnatlon 
























RESIDUAL STRENGTH SERVICE DATA 

SERVICE LIFE 
INSPECTION INTERVAL 
REPAIR 

Figure 1. - CODSTRAN computer code schematic. 



Figure 2. - CODS FRAN flow chart: computational simulation of progressive fracture. 




Igure 4. - Composite laminate finite element model (approximately 
540 elements, 450 nodes, 6/0 DOf ), 
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